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ABSTRACT 

Tiny grains such as polycyclic aromatic hydrocarbons (PAHs) have been thought to dramatically 
reduce the coupling between gas and magnetic fields in weakly ionized gas such as in protoplanetary 
disks (PPDs) because they provide tremendous surface area to recombine free electrons. The presence 
of tiny grains in PPDs thus raises the question of whether the magnetorotational instability (MRI) 
is able to drive rapid accretion to be consistent with observations. Charged tiny grains have similar 
conduction properties as ions, whose presence leads to qualitatively new behaviors in the conductivity 
tensor, characterized by n/rie > 1, where Ue and n denote the number densities of free electrons 
and all other charged species respectively. In particular, Ohmic conductivity becomes dominated by 
charged grains rather than electrons when n/n^ exceeds about 10^, and Hall and ambipolar diffusion 
(AD) coefficients are reduced by a factor of {n/rie)'^ in the AD dominated regime relative to that in 
the Ohmic regime. Applying the methodology of Bai (2011), we find that in PPDs, when PAHs are 
sufficiently abundant (> 10~^ per H2 molecule), there exists a transition radius rtrans of about 10 — 20 
AU, beyond which the MRI active layer extends to the disk midplane. At r < rtrans, the optimistically 
predicted MRI-driven accretion rate M is one to two orders of magnitude smaller than that in the 
grain-free case, which is too small compared with the observed rates, but is in general no smaller 
than the predicted M with solar-abundance 0.1/^m grains. At r > rtians, we find that remarkably, 
the predicted M exceeds the grain-free case due to a net reduction of AD by charged tiny grains, and 
reaches a few times 10~*AfQ yr~^. This is sufficient to account for the observed M in transitional 
disks. Larger grains (> 0.1//m) are too massive to reach such high abundance as tiny grains and to 
facilitate the accretion process. 

Subject headings: accretion, accretion disks — instabilities — magnetohydrodynamics (MHD) — pro- 
toplanetary disks — turbulence 



1. INTRODUCTION 

Weakly ionized gas is not perfectly conducting. In 
the presence of magnetic field (B), the Ohm's law is 
anisotropic, giving rise to Ohmic resistivity. Hall effect, 
and ambipolar diffusion (AD) as the three non-ideal mag- 
netohydrodynamics (MHD) effects. More formally, as- 
suming the inertia in the ionized species is negligible, 
the dynamical equations for the neutrals are the same 
as ideal MHD except for the induction equation, which 
reads 



reads (jSalmeron fc Wardlel[2003[) : 
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(1) 
where v is the velocity of the neutrals, J ~ (c/47r) V x B 

is the current density, B denotes unit vector along B, 
subscript "^" denotes the vector component that is per- 
pendicular to B, rjo, TlH and tja arc the Ohmic, Hall and 
the ambipolar diffusivities, which depend on the num- 
ber density of the charged species, and the latter two 
also depend on B. When ions and electrons are the 
only charged species in the gas (which can be relaxed 
to when the electrons are the main negative charge car- 
riers), a simple expression for the diffusion coefficients 
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where rje is the resistivity due to the electrons, rig is the 
electron number density, p and pi are the mass density of 
the neutrals and the ions, 7^ is defined after equation ^. 
In general, Ohmic effect dominates in dense regions with 
weak magnetic field, AD dominates in tenuous regions 
with strong magnetic field, while the Hall regime lies in 
between. 

Grains, which typically contribute to about < 1% of 
mass in the gas, can possess both positive and nega- 
tive charges. Grains are well known for significantly en- 
hanc ing the electron recombination rate (jDraine fc SutinI 
Il987| ). which reduces the ionization fraction {ne/riH, 
where nn is the number density of hydrogen atom) , and 
leads to greatly increased Ohmic resistivity. However, 
in sufficiently weakly ionized gas, the abundance of tiny 
grains (e.g., the polycyclic aromatic hydrocarbon, PAH) 
may well exceed the electron abundance without con- 
tributing much to the mass budget. In this paper, we 
generalize equation ([2|) to include positive and negatively 
charged tiny grains and show that for sufficiently weak 
ionization, tiny grains can dramatically reduce r]H and 
rjA as the main charge carrier switches to the grains in- 
stead of ions and electrons. 

Our result is directly applicable to protoplanetary disks 
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(PPDs) , where the ionization fraction rcsuhing from ma- 
jor ionization sources such as X-rays and cosmic rays is 
generally orders of magnitude below unity, and it de- 
creases from surface to the mid plane because o f the at- 
tenuation of ionizing particles (jGammid Il996| ) . Non- 
ideal MHD effects are closely relevant to the rnagne - 
torotational instability (MRI, iBalbus fc HawlevI [l99ll ). 
whose linear dispersion relation as well as t he non- 
linea r saturation pr operties are strongly affected (jWardld 
Il999t lBalbuj[2009l) . and whether the MRI is responsi- 
ble for driving rapid accretion with M ~ 10~^^^A'jfQ 
yj.-i jjj^ PPDs (jHartmann et al.lll998[ ) has been a long- 
standing problem. Most studies h ave focused on the 
role of the Ohmic res istivity (e.g., I Turner et al.l 120071 : 
iBai fc GoodmanI |2009|) ■ while it was not until recently 
have AD been taken int o account to estimate the ef- 
fectiveness of the MRI (IChiang fc Murrav-CTavl 120071: 
iPerez-Becker fc Chiang||2011a^ . Combining the most re- 
cent results from numer ic al simulati ons of the MRI with 
AD (jBaifc Stond 120111 ) ■ iBal (j2011h showed that while 
the MRI can always operate in PPDs, AD can be a 
main limiting factor for the MRI-driven accretion (be- 
sides Ohmic resistivity) . The reduction of AD coefficient 
by tiny grains (e.g., PAHs), while counterintuitive, helps 
enhance the accretion rate, as we will demonstrate in this 
paper. 

Although observational data show strong evidence 
of gra in growth to micron size or larger in PPD s 
(e.g., iD'Alessio etldl l2001l : Ivan Boekel et al.l I2003D . 
PAH emission has also bee n detected in major- 
ity o f Herbig Ae/Be stars (jAcke fc van den Anckerl 
|2004|). as well as a small fraction o f T-Tauri stars 
teecrs ct al.' '2006'; 'Olivcira et alll2010( ). As argued in 
[Pcrez-Bccker & Chiang (20 11al ). PAHs may be equally 
abundant in T-Tauri disks but they fluoresce less lumi- 
nously due to fainter ultraviolet radiation field of their 
host stars. The existence of PAHs also suggests a contin- 
uous size distribution of grains to the smallest end of a 
few A as a result of grain coagulation and fragmentation. 
Throughout this paper, we use the phrases "tiny grain" 
and PAH interchangeably, which refer to grain with size 
a < 0.01/im. Tiny grains may dominate larger grains in 
abundance while contribute a negligible fraction of the 
total grain mass. Taking the grain density to be 3 g 
cm~^ and the gas mean molecular weight /i„ to be 2.34 
atomic mass (as appropriate for PPDs), the relation be- 
tween grain mass fraction (/) and grain abundance per 
H2 molecule [x) reads 

^ ~0.32f^V:..:^V, (3) 



0.01 VlO-V VO-OlAim, 

where a is grain size. The abundance of grains with 
a < 0.01/im may easily reach large abundance of 10~^ or 
higher, while the abundance of a > 0.1 /xm grains would 
be at most about 10~^^. The significance of this abun- 
dance cut {x = 10~^) will be addressed in this paper. 

In Section [2] we describe a generalized model for non- 
ideal MHD diffusion coefficients with the inclusion of 
charged grains. The model is applied to interpret the 
the results in Section [3l where we study the role of tiny 
grains in PPDs following the methodology of IBai, (20 11.1 . 
Summary and discussion follow in Section |4l 



2. NON-IDEAL MHD EFFECTS WITH GRAINS 

The Ohm's Law derives from the motion of charged 
particles, which is characterized by the Hall parameter, 
the ratio betwe en the gyrofre quency and the momentum 
exchange rate (|Wardldl2007[ ). For species j with mass 



and charge Zje, the Hall parameter reads 
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where 7^ = {<7v)j/{fin+'nij) with (cru)j being the rate co- 
efficient for momentum transfer between charged species 
j with the neutrals and fin is the mean molecular weight 
of the neutrals. Charged species j is strongly coupled to 
the neutrals if |/3j| <ti 1, and is strongly tied to magnetic 
fields when |/3j| > 1. 

In weakly ionized gas with a number of different 
charged species, the general expressions for the Ohmic, 
Hall and ambipolar diff usion coefficients arc given by 
(|Wardldl200"7t iBallMll) 
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where a± = vOff+^p ^^"^ the Ohmic, Hall and Peder- 
sen conductivities are 
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respectively, where the summation goes over all charged 
species. Note that the Hall conductivity depends on the 
sign of Zj, while ap and ap depend only on \Zj\. 

The simple expressions for the magnetic diffusion co- 
efficients (m can be obtained by assuming electrons and 
ions arc the only two charged species. Below we gener- 
alize it to include tiny grains. Besides the reason men- 
tioned in the introduction, we consider tiny grains be- 
cause they are unlikely t o possess multiple charges due t o 
higher potential barrier ([Perez-Becker fc Chianeir2011al) . 
which simplifies the algebra considerably. In addition, 
the mean free path for tiny grains is sufficiently small so 
that they can be treated as fluid. 

The momentum transfer rate coefficients for electrons, 
ions and grains can be found in equations (14) - (16) of 
[Bai (2011). Converting to the Hall parameter, we have 
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for ions and electrons respectively, where Bp is the mag- 
netic field strength measured in Gauss, nis = n/f/lO^^ 
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cni^ , and the /i„ is taken to be 2.34 atomic mass ap- 
propriate for PPDs. For singly charged smaU grains, 
the resulting (3 is very close to /3i, because the reduced 
gyro- frequency due to larger mass is compensated by the 
reduced momentum transfer rate due to larger inertia. 
The electron Hall parameter /3e is much larger than Pi 
because electrons are the most mobile species. There- 
fore, in terms of conductivity, we only have two groups 
of charges: electrons with Hall parameter /?£, and ions / 
charged grains with Hall parameter /3i. 

Let n = rii + n^^ + n~ be the total number density 
of the second group. Since grains can carry negative 
charge, one has fi > n^. Generalization from n = Ug^ 
(grain-free case, which leads to equation ([2])) to fi > n^ 
involves more complicated algebra, as we describe below. 
For brevity, we ignore the pre-factor eriec/ B in conduc- 
tivities and cB/ATTerif, in diffusivities, which make them 
dimensionless. 

The three components of the conductivity tensor read 
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The parameter 9 is independent of magnetic field, and 
measure the ratio of grain conductivity to electron con- 
ductivity. Note that = in the grain-free case, and the 
Hall conductivity is independent of 9. The perpendicular 
conductivity is 



where 



^^ = 



fiO)^ 



We + Pi)^ 



(1 + /3D 



[{l + 9f + {P, 



^Y 



(10) 



(11) 



Note that f{9) = 1 in the grain- free case {9 = {)). It is 
straightforward to obtain the Ohmic, Hall and ambipolar 
diffusion coefficients, after some algebra 
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where the approximate formulae are obtained by noting 
that Pe^ Pi, with error on the order of Pi/Pe ^ 10"'^. 



We see that the grains affect the magnetic diffu- 
sion coefficients mainly via two factors: (1 + 6) and 
(Pi + 9Pf.)/ Pi ~ n/ue, both of which are independent 
of magnetic field strength. The former describes the ra- 
tio of total conductivity to electron conductivity, while 
the latter describes the number density ratio of ions / 
grains to electrons, and is much more substantial than 
the factor (1-1-6') since Pe':^ Pi- 

For 9 ^ Pi/Pe ^ 10~^, we generally have n^ rii k, rie, 
f{9) « 1, and the classical results ^ are recovered. In 
particular, Ohmic resistivity dominates when /3e < 1, 
where both ions and electrons are coupled to the neutrals; 
Hall regime corresponds to Pi < 1 < Pe, where electrons 
are coupled to the magnetic field while the ions are not; 
AD regime corresponds to Pi > 1 where both electrons 
and ions are coupled to the magnetic field. This is the 
classical interpretation of the non-ideal MHD effects. 

Qualitatively new behaviors appear when 9 > Pi/Pe-, 
where grains play a dominant role in the abundance of 
charged particles (nE. ^ ng). In the next two paragraphs, 
we discuss the asymptotic behaviors of the diffusion co- 
efficients, where comparison is made with the grain- free 
diffusion coefficients ([2]) at fixed n^- In reality, tiny grains 
strongly reduces Tig relative to the grain-free case, which 
will be discussed at the end of this section. 

The Ohmic resistivity is least affected by grains, which 
is simply reduced by a factor of (1 + 6') relative to the 
electron resistivity regardless of magnetic field strength. 
Therefore, grain conductivity becomes important when 



> 1 or n > 10-^ 



For Hall and AD coefficients, we 



consider two separate limits. When /3i ^ /3e ^ 1 (Ohmic 
regime), we have 
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where we have factored out the results into the grain-free 
expression (left) multiplied by a correction factor (right). 
We see that the Hall diffusivity is only moderately af- 
fected by grains, being reduced by a factor of (1 -I- 9)'^ 
at fixed TT-e. The AD, on the other hand, is enhanced 
by a factor of n/rie, in addition to a moderate reduction 
by (1 -t- 9Y . Nevertheless, changes in Hall and AD co- 
efficients do not play a significant role here since Ohmic 
resistivity is still the dominant effect. 

In the opposite limit 1 <^ Pi <^ Pe (AD regime), we 
have 
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We see that in this limit and at fixed n^, the Hall effect 
is reduced by a factor of {fi/rieY, while AD is reduced 
by a factor of in/ne). The reduction of AD is easily un- 
derstood. Without grains, the neutral gas is coupled to 
the magnetic field through ion-neutral collisions, hence 
rjA oc \/ni = l/rig. Electrons play a negligible role be- 
cause of its small inertia. Charged grains (no matter 
positive or negative) play exactly the same role as ions. 
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Fig. 1. — Dimonsionlcss Ohmic (black), Hall (blue) and ambipo- 
lar (red) diffusion coefficients as a function of the electron Hall 
parameter /3e = 1000/3^. We have considered two different values 
of the charged grain abundance parameter: d = 10"'^ (n = 2n^, 
dashed) and 9 = 0.25 (n = 251ne, solid). 

hence we have tja oc \/n in the presence of grains. We 
also note that the Hall effect vanishes if positive and 
negative charge carriers have the same mass, which is 
consistent with our result as rie/n — > 0. 

In Figure [T] we show two sample calculations of the 
magnetic diffusion coefficients for = 10""^ (dashed) and 
= 0.25 (solid) respectively, which corresponds to sit- 
uations where charged grain abundance is comparable 
to, and greatly exceed the electron abundance. The two 
asymptotic regimes derived above are clearly seen, with 
a transition region in between. The 9 = 10"'^ results 
are close to the grain-free case, where all three curves 
are close to straight lines except weak transitions as Pi 
passes 1. For relatively large 6*, the transition region ex- 
tends from /3e Ri 1 to /3i « 1. Because of the suppression 
of the Hall diffusion beyond /3e ~ 1 and enhancement of 
AD before /3e ~ 1, AD becomes the dominant effect even 
when Pi <^1 [jSe <Si 1000), and the Hall regime gradually 
diminishes as 6 increases. These results nicely explain the 
observ ed magnetic diffusivity pattern in Figure 2 of lBail 
(|2011f ). as aided with the chemical abundances shown in 
Figure 1 of the same paper. 

Now let us take into account the chemistry in the 
gas. Assuming fixed gas density, temperature, and ion- 
ization rate, the equilibrium electron number density in 
the presence of small grains n^i is much smaller than 
that in the grain-free case We p ()Bai fc GoodmanI 120091 : 
iPerez-Becker fc Chiandl2011aD . Therefore, Hall and AD 
coefficients are made larger by a factor of neo/j^ei- How- 
ever, this factor is compensated by another reduction 
factor of {nei/nY' for the Hall effect and of Uei/n for 
AD in the AD dominated regime. Here we focus on AD. 
In order to have a net reduction of AD coefficient when 
/3i » 1 (AD regime), it is required that 

n > neo ■ (17) 

That is, the charged grain abundance has to exceed the 
grain-free electron abundance. We will see that this con- 
dition can be met in PPDs. Similarly, for Ohmic resis- 
tivity to be reduced relative to the grain-free case, one 



requires n > lO^neo- This condition, however, is very un- 
reahstic given the disk chemistry (see Section [3T|l . thus 
Ohmic resistivity always increases in the presence of tiny 
grains. Nevertheless, if condition p7|) is met (which is 
generally the case when PAHs are sufficiently abundant, 
see Section r3.1|) . Ohmic resistivity is enhanced by at most 
a factor of about lO'^ relative to the grain- free case. 

3. APPLICATION TO THE MAGNETOROTATIONAL 
INSTABILITY IN PPDS 

To illustrate the significance of the theoretical results 
in the previous section, we study the role of tiny grains 
in PPDs. We refer to Section 3 of JBa] (|2011f ) for all de- 
tails about the methodology and calculation procedures, 
which we follow exactly except that here we consider 
a = 1 nm sized grains (as a proxy for PAHs) instead of 
larger gra his^ . The PAH abu ndance in T-Tauri stars esti- 
mated bv lGeers et all (12006') is a bout xpah ~ 10~^-10~J 
per H2 molecule, while Perez-B ecker fc ChiansJ ()2011al ) 
argued for smaller abundance due to dust settling. More- 
over, the distribution of PAHs in PPDs ha s been found 
to be spatially variable (jGeers et al.ll2007t ). For these 
reasons, we adopt ccpah = 10~® as fiducial, while we also 
consider a;pAH down to lO""*^""^ and up to 10~^. Wc fix the 
X-ray luminosity to be lO'^*^ erg s^^, X-ray temperature 
to be 5 keV, and include cosmic-ray ionization of 10^^'' 
s~^ with penetration depth of 96 g c m~^. We use the 
minimum-mass solar nebular (MMSN, lHavashilll98l[ ) as 
the fiducial disk model, while we als o cons ider the much 
more massive disk model bv iDeschI (l2007l ). A complex 
chemical network of more than 4000 reactions are evolved 
at different heights and radii of the disk for 1 Myr to 
(quasi-) chemical equilibrium, from which abundance of 
charged species are extracted to evaluate the non-ideal 
MHD diffusion coefficients. We will mainly focus on the 
outer disk (> 10 AU) where AD rather than Ohmic re- 
sistivity plays the dominant role the gas dynamics and 
the effect discussed in Section [2] is most relevant. 

3.1. Abundance of Charged Species 

In Figure [2] wc show the electron, ion and grain abun- 
dances at 40 AU in the MMSN model as a function of disk 
height (z), normalized to the disk scale height {H = c^/O 
where Cg is the sound speed, fi is the disk angular fre- 
quency). We see that electrons and ions are the domi- 
nant charged species (i.e., n k n^) above ^ 3H, while 
below it the abundances of ions and charged grains ex- 
ceed electron abundance, and n/ue reaches about 1000 
at the disk midplane. The large value of n/ue implies 
strong supprcsion of AD and the Hall effect according to 
equation p^ . To see its significance, we also show the 
electron (thus ion) abundance in the grain-free calcula- 
tion as the red dash-dotted line in the Figure. Clearly, 
we see that Ueo < n for z < 2H . According to equation 
p7)) . this means that when tiny grains are present, the 
AD coefficient is even smaller than that in the grain-free 
case! Seemingly counterintuitive, one can qualitatively 
understand this result as follows. 

^ We use binding energy of D = 3 eV in the calculatio n of 
the el ectron-grain sticking coefficient rather than 1 eV used in lBail 
1I20IIII 
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Fig. 2. — Abundance of charged species and grains (PAHs) as 
a function of disk height in our fiducial model calculation (MMSN 
disk with xpAH = 10~*) at 40 AU, plotted as solid lines with 
labels. For comparison, red dash-dotted line shows the electron 
(and ion) abundance in the grain-free calculation. The ionization 
rate is plotted (black dashed) for reference. 



Near the disk midplane, we see that rigr 3> n. 



± > 

gr 



Tlj > 



Ue- In this regime, the electrons produced from H2 ion- 
ization is quickly swallowed by the grains. Similarly, the 
ions exchange charge with neutral grains to produce pos- 
itively charged grains. Therefore, ionization effectively 
takes place on the grains: 2 gr — )■ gr"*" -I- gr~, with the 
same ionization rate <^eff- The dominant recombination 
channel is simply its inverse re action, with recombina- 
tion ra te given by equation (3) of lUmebavashi fc Nakand 
(|1990f l. For recombination between two equal sized tiny 
grains {e'^/akT 3> 1), it reduces to 



48e^ 
PciakT 
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where pd = 3 g cm^^ is the grain mass density, k is the 
Boltzmann constant, T is the temperature. In the second 
equation ai denotes grain size normalized to 1 nm, and 
T'loo = 2^/100 K. As long as charged grain recombination 
is the dominant recombination process, the abundance 
of charged grains can be approximately given by 
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where Ccff, -17 is the ionization rate normalized to 
10"^^ s~^, riH, 10 is the number density of the hydrogen 
atoms normalized to 10^*^ cm^"^. Plugging in the num- 
bers relevant to Figure[2]at midplane, we find n±./nH ~ 
1.1 X 10^^". Our chemistry calculation gives 6.5 x 10"^^ 
for the averaged abundance of charged grains, which is 
slightly smaller due to small contributions from other re- 
combination channels, but is within a factor of 2 from 



the analytical estimate. Our chemistry calculations fur- 
ther reveal that n±./nH increases weakly with total grain 
abundance rigr (or xpah) and approaches the asymptotic 
value (fT9|). 

In the grain-free case, the electron abundance rieo is 
determined by the balance between ionization and mul- 
tiple recombination channels, dominated by dissociative 
recombinations. Typical electron-ion dissociative recom- 
bination rate coefficients are on the order of 10~^ s~^ 
cm^^ at lOOK, which is a factor of several higher than the 
grain recombination coefficient (|18|) . The higher recom- 
bination rate leads to smaller ionization level than our es- 
timate (|19p . which explains why n^o < fi in the presence 
of abundant tiny grains. We note that the value of rieo 
depends on the choice of chem ical reaction network. Sim- 
ple re action network (such as lOppenheimer fc Dalgarng 
Il974f l generally produces larger neo mainly because of 
the lack of recombination channels. In our complex 
(and presumably more realistic) network, we find that 
the dominant recombination process is due to NII4 and 
CII3CNH+ in this particular case, giving neo/n^ to be 
about 3.2 X 10^^^. As a result, we obtain fi/yieo ~ 4 in 
the midplane, which leads to a substantial net reduction 
of the AD coefficient. 



3.2. Active Layer and Accretion Rate 

The fact that n > ?T,eo implies that tiny grains may 
facilitate the MRI by suppressing AD. To see this more 
explicitly, we calculate the non-ideal MHD diffus ion coef- 
ficients and apply the criteria (20) in iBail (|2011f ) to iden- 
tify the MRI-activc regions in our adopted PPD model. 
Briefly, for the MRI to operate, the Ohmic Elsasser num- 



ber A = 



jfi has to be greater than unity, where 



VA = Bj \JAiip is the Alfven velocity. Moreover, the 
magnetic field has to be weaker than some critical value, 
which is set by /3 > /3min(Am). Here plasma /3 is the ra- 
tio of gas to magnetic pressure (not to be confused with 
the Hall parameter), Am = v'j/riA ^ is the AD Elsass er 
number, and f3,^in{Am) is given by (iBai fc Stonell2011[ ) 
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(20) 



which increases with decreasing Am. For the number 
density profile given in Figure [2l the result is shown in 
Figure [H 

We see that the inclusion of tiny grains strongly in- 
creases the Ohmic resistivity, as expected. At disk mid- 
plane, magnetic field has to be about 10 times stronger 
than the grain- free case in order for A to be above 1 (i.e., 
rjo is about 100 times larger^). However, the ionization 
fraction at disk midplane is still large enough such that 
in both cases, the required field strength for A > 1 is 
still far from equipartition (i.e., /3 ^ 1). This fact makes 
Ohmic resistivity essentially irrelevant in determining the 
vertical extent of the active layer and the accretion rate 
(since the midplane is already active), while the decisive 
role is played by AD. 

^ Note that since n fa lO^Ue as read from Figure [2] electrons 
and other charged species contribute roughly equally to the Ohmic 
conductivity. 
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Fig. 3. — MRI permitted region (sliadcd) in the MMSN model at 40 AU from the grain-free calculation (left) and the calculation 
including PAHs with xpAH = 10~® (right). Bold black solid line marks the boundary given by the Ohmic Elsasser number criterion A = 1. 
Thin black solid lines show contours of constant AD Elsasser number Am, labeled by log j^q (Am). Bold and thin blue lines mark plasma 
/3 = 1 and /3 = 100 respectively. The red dashed lines indicate the required field strength in the MRI permitted region in order for driving 
accretion rate of lO^'^, 10^* and IO^^A/q yr~^. 

In Figure [31 contours of constant Am are straight lines 
in the grain free case, as Am is independent of the mag- 
netic field. With the inclusion of tiny grains, we see that 
near the midplanc region. Am increases by 3 orders of 
magnitude as magnetic field increases. This corresponds 
to the transition region in Figure [T] with Pi < 1 < Pe- 
Eventually, at sufficiently strong magnetic field {Pi > 1), 
Am, becomes independent of magnetic field again. Since 
h > TT-eo, Am is larger in the presence of PAHs than in the 
grain- free case. Therefore, according to (f20| . stronger 
magnetic field is permitted near the disk midplane thanks 
to the PAHs. According to the strong linear correlation 
between the Shakur a-Sunyaev a pararti eter and the mag- 
netic field strength (jBai fc Stonell2011l ). an approximate 
magnetic field strength can be derived if MRI operates 
in PPDs (jBaill201lD 



B « 2dMn/H ' 



LOM^/VX"/^ G 



(21) 



which assumes the thickness of the active layer to be H 
on each side, and the second equation assumes one solar 
mass protostar with MMSN temperature profile. Here 
M_8 is accretion rate measured in IO^^A/q yr"-'^, tau 
is disk radius measured in AU. The required magnetic 
field strengths for typical observed accretion rates are 
also illustrated in the Figure. Since stronger field leads 
to faster accretion, it becomes clear that PAHs are able 
to make PPDs accrete more rapidly than in the grain-free 
case. 

We note that the above explanation for the enhance- 
ment of accretion by tiny grains no longer completely 
holds in the inner region of PPDs (< 10 AU), where the 
disk surface density is far above the X-ray penetration 
depth and accretion is layered. In such situations, the 
lower boundary of the active layer is determined by both 
Ohmic resistivity and AD for the two Elsasser number 
criteria to be satisfied (e.g., see the left panel of Figure 
[S] in the Appendix) . Since Ohmic resistivity is larger in 
the presence of grains (see the end of Section [2]), accre- 



tion is thus less efficient than the grain-free case. On the 
other hand, when PAHs are sufficiently abundant, the 
increase in Ohmic resistivity is bounded (< 10'^) since 
conductivity is dominated by the charged PAHs rather 
than electrons, and the reduction in accretion rate is in 
general no more than the case with 0.1/im grains (see the 
next subsection). 

3.3. Parameter Study 

A more qua ntitative es timate of the disk accretion rate 
is discussed in lBail ()2011f ) (see their Section 4.3). The pre- 
dicted accretion rate M is based on the hypothesis that 
MRI tends to arrange the field configuration to maxi- 
mize the rate of angular momentum transport, and is 
optimistic. To demonstrate the role of tiny grains more 
extensively, we perform our calculations for both the 
MMSN model and the more massive Desch's disk model 
at a wide range of disk radii from 1 AU to 100 AU, and 
consider PAH abundances of xpah = 10~^, 10~® (fidu- 
cial) and 10~^. In addition, we conduct one calculation 
with cosmic-ray ionization turned off at xpah = 10~^ for 
both disk models. In Figure SI we show the predicted M 
as a function of disk radius for all our calculations. This 
Figure is to be compared with Figure 4 of lBai (2011i) , who 
considered the cases with 0.1 ^m and 1/xm grains. We see 
that the inclusion of tiny grains makes the predicted M 
comparable to that with solar abundance 0.1/im grains 
at 1 AU, which is one to two orders magnitude smaller 
than the grain-free case. However, at larger disk radii, 
tiny grains promote accretion compared with sub- micron 
grains, with predicted rate that even exceeds the grain- 
free case. Below we discuss the effect of PAHs in more 
detail. 

We find that whenever the active layer extends to the 
disk midplane (the gray area touches the vertical axis 
in Figure [3l which makes the Ohmic resistivity irrele- 
vant), the predicted M in the presence of PAHs exceeds 
the grain-free case, and the reason is simply due to the 
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Fig. 4. — Optimistically predicted accretion rate as a function of disk radius for the MMAN disk (left) and the Dcsch's disk (right) models. 
In each panel, we show results from grain-free calculations (black circles) and calculations with various PAH abundances as indicated in 
the legend. For xpah = 10~*, we also show the results from calculations without cosmic-ray ionization in the black dash-dotted line. The 
cross symbols in each panel denote results from the same set of calculations at selected radii but with xpAH = 10~^^. 

reduction of AD coefficients as discussed in the previous 
two subsections. This occurs at disk radii r > rtmns ~ 15 
AU for xpAH = 10""^ and 10~^, and for r > rtrans ~ 40 
AU for ccpAH = 10^^. Inside the transition radius rtj-ansj 
accretion is layered and the predicted M is weU below 
that in the grain-free case due to Ohmic resistivity as 
discussed at the end of the previous subsection. The in- 
crease in predicted M near rtians is very sharp, which is 
closely related to dependence of Am on magnetic field 
strength shown in Figure [1] and is discussed in more de- 
tail in the Appendix. 

An interesting and counterintuitive fact is that higher 
PAH abundance leads to faster accretion at all disk radii. 
This fact can be understood as follows. First of all, 
near the base of the active layer, we find that charged 
grain abundance is so much higher than the electron 
abundance that the Ohmic resistivity is determined by 
charged grains rather than electrons (i.e., 6 > 1). There- 
fore, even the abundance of free electrons rapidly de- 
creases with increasing ccpah, the Ohmic resistivity ap- 
proximately stays at the same level (while still much 
larger than the grain- free resistivity). Secondly, as we 
mentioned in the discussion after equation (|19p , the grain 



ionization level n^^/nH increases weakly with the to- 
tal grain abundance xpah- This fact reduces both the 
Ohmic resistivity and the AD coefficient, which leads to 
faster accretion. Only at r < 1 AU, and when xpah is 
small (10~^ or less), does the Ohmic resistivity at the 
base of the active layer determined by electrons, which 
produces higher accretion rate than larger zpah cases. 

We see that in the presence of PAHs, the predicted M 
increases with disk radius for 



< 



rtr 



where accretion 
> 



is layered, and falls off with radius for r > rtrans as the 
disk midplane is activated. Also, larger disk surface den- 
sity generally leads to smaller M in the inner disk when 
accretion is layered (although not by much) , while when 
the disk midplane becomes active in the outer disk, the 
Desch's disk model gives higher M. These features are 



all consistent with the results presented in iBail ()2011|) , 
where explanations are offered (see Section 5.1). 

In addition, we see in the dash-dotted line of Figure 
m that when cosmic-ray ionization is turned off (since 
it may be shielded by protostellar winds), the predicted 
M is almost indistinguishable with our fiducial case at 
the inner 10 AU. The predicted M still exceeds that in 
the grain-free case in the outer disk, with the transition 
radius rtrans slightly larger than that in the fiducial cal- 
culation, and the increase in M near rtrans is less sharp. 
Beyond rtrans, the accretion rate reduction is generally 
within a factor of 2 of the fiduci al result. These results 
are all consistent with iBail (|2011[ ) , where it is shown and 
discussed that the deeper-penetrating cosmic-ray ioniza- 
tion does not substantially enhance the accretion rate 
compared with the X-ray ionization (see Section 5.2). 

Finally, we emphasize that the existence of rtrans and 
the fact that M increases with xpah holds only when 
PAHs are sufficiently abundant (xpah ^ 10~^). For 
comparison, we also show in Figure |4] with cross sym- 
bols the predicted accretion rate at selected disk radii 
with xpAH = 10"^"'^. At small radius of ~ 1 AU, the 
predicted M is well above the cases with higher PAH 
abundances. This is because the free electrons is much 
more abundant so that electron conductivity dominates 
grain conductivity, and gives smaller resistivity. The disk 
midplane becomes active at r > 15 AU, but we do not 
find the sharp jump in the predicted M as in the case 
with higher PAH abundances, and the predicted accre- 
tion rate is also smaller by a factor of up to 100. This is 
because xpah is so small that the grain number density is 
well below rie (so h/rie ~ 1) and the effects we discussed 
before no longer apply. 

The role of PAHs in PPDs has be en recently stud- 
ied by iPerez-Becker fc Chiand (|2011aD . Besides the fact 
that they used a more stringent criterion for AD to 
suppress the MRI, the Am values they obtained were 
largely underestimated because they adopted the sim- 
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Fig. 5. — Optimistically predicted accretion rate as a function 
of disk radius for the MMSN disk with two populations of grains; 
0.1/^m grains with fixed mass fraction of 1% and PAHs with abun- 
dance xpAH = 10~^,10~* and 10"'^. Comparison is made to 
the grain-free case (black solid) and case with single population 
of 0.1/im grains (black dash-dotted). 

plc foriTLulac ([2]) without taking into account the con- 
tribution from grains. This has led them to conclude 
that the presence of PAHs makes the MRI too inefficient 
to drive rapid accretion at all disk radii, although as 
we have shown that PAHs in fact promote accretion in 
the outer disk, and at all disk radii the predicted ac- 
cretion rate generally increases with PAH abundance. 
Our findi ngs are subsequently confirmed in their follow- 
up work (jPerez-Becker fc Chiangir2011b[ ) with corrected 
grain conductivities. 

3.4. Effect of Grain Size Distribution 

Our study so far have considered only single-species of 
grains. In reality, tiny grains coexist with larger grains 
with the latter dominates the grain mass (but not nec- 
essarily abundance). To justify our findings in a more 
realistic context, we conduct another set of calculations 
with two populations of grains. Besides the PAH popula- 
tion with a;pAH = 10^^, 10~^ and 10~^ that contributes a 
negligible amount of mass, we add the second population 
of 0.1/im grains with fixed mass fraction of 1%. In Figure 
[31 we show the predicted accretion rate as a function of 
disk radius in the same way as Figure |4l For compari- 
son, we also plot the predicted M in the grain- free case 
(black solid) and the pure O.l^m grain case (black dash- 
dotted) taken from Bai (20li|). We see that at fixed PAH 
abundance with xpah ^ 10~^, the predicted M with two 
grain populations is very close to that in our previous cal- 
culations with a single PAH population, and it behaves 
as if the 0.1/im grains do not exist. In particular, the 
sharp enhancement of accretion beyond the transition 
radius is clearly present, and the predicted M increases 
with .TpAH- At xpAH = 10~^, the situation is marginal: 
with the addition of 0.1/im grains, the signature of sharp 
transition is less obvious, which occurs at rtrans ^ 60 



AU, and beyond rtrans the predicted M is still less than 
that in the grain-free case. In fact, the behavior of the 
two-population calculation at xpah = 10~^ roughly lies 
in between the two single-population limits. 

Our results ar e consi stent with chemistry studies by 
iBai fc GoodmanI (j2009t ) which indicate that different 
grain populations tend to behave independently. There- 
fore, the outcome of a size distribution of grains is mainly 
determined by the grain population that d omina tes the 
recombination process. iBai fc GoodmanI (I2009D found 
that the controlling parameter of the recombination lies 
in between the total surface area and the grain abun- 
dance weighted by linear size. As an approximation we 
weigh the grain abundance by a?/'^ and find that PAHs 
and 0.1/im sized grains contribute roughly equally to the 
recombination process when a;pAH ^ 3 x 10^^, consistent 
with our findings discussed in the previous paragraph. 

In sum, the effect of tiny grains studied in this paper 
is approximately independent of the presence of larger 
grains, and we quote xpah = 10~^ as the critical tiny 
grain abundance above which enhancement of accretion 
is possible in the outer regions of PPDs. 

4. SUMMARY AND DISCUSSION 

Tiny grains are important charge carriers in weakly 
ionized gas and strongly affect the non-ideal MHD ef- 
fects when their abundance is above the ionization frac- 
tion. We generalize the commonly used simple expres- 
sions for non-ideal MHD diffusion coefficients ^ to in- 
corporate the effect of charged tiny grains. We show 
that at sufficiently small ionization level, tiny grains be- 
come the dominant charge carriers, and the non-ideal 
MHD diffusion coefficients behave very differently from 
the grain-free case. The Ohmic conductivity is domi- 
nated by charged grains rather than electrons when n 
exceeds about 10"^ne. In the AD regime (strong magnetic 
field). Hall and AD coefficients are strongly reduced by 
a factor of about (n/rieY relative to those in the Ohmic 
regime (weak magnetic field), and for sufficiently large 
n/rie^ Hall dominated regime diminishes. 

We study the role of tiny grains on the MRI driven 
accretion in PPDs, and find that novel behaviors oc- 
cur when the tiny grains are sufficiently abundant with 
a^PAH ^ 10~^, regardless of whether larger grains are 
present or not. At the inner disk where accretion is lay- 
ered, the predicted accretion rate in the presence of tiny 
grains is one to two orders of magnitude less than the 
grain-free case due to increased Ohmic resistivity, but 
is similar to or higher than that with solar-abundance 
0.1/im grains. A sharp increase in the predicted M occurs 
at the transition radius rtrans ~ 15 AU (in the fiducial 
model) where the disk midplane becomes active, making 
Ohmic resistivity irrelevant to the accretion rate. Quite 
unexpectedly, we find that at r > rtrans, tiny grains make 
accretion even more rapid than the grain-free case. More- 
over, our predicted accretion rate increases with PAH 
abundance. These results are due to that at large PAH 
abundance, ionization-recombination balance makes n 
orders of magnitude larger than rig, and even exceeds the 
grain-free electron density n^Q at disk midplane. These 
facts prevent Ohmic resistivity from rapidly increasing 
as xpAH increases, reduce the dissipation by AD, thus 
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Fig. 6.— Same as Figure |3] but for MRI permitted region (shaded) in the MMSN model at 10 AU (left) and 16 AU (right) from 
calculations with PAHs at xpAH = 10^**. 



facilitate the active layer to extend deeper into the disk 
midplanc and permit stronger MRI turbulence. Our re- 
sults highlight the importance of evaluating the full con- 
ductivity tensor in the calculation of the non-ideal MHD 
diffusion coefficients rather than using the simple grain- 
free formulae ([2]). 

We emphasize that the effects studied in this paper 
mainly apply to tiny grains (< 0.01/im), and when they 
are sufficiently abundant (xpah ^ 10~^). For grains 
larger than O.lfim., their abundance is at most 3 x 10~^^ 
per H2 molecule (for having 1% of mass), which is or- 
ders of magnitude smaller than 10~^. The reduction of 
AD coefficient also exist for these relatively larger grains, 
as can be seen in the bottom panels of Figure 3 in iBail 
(|2011h . but its effect is much more limited than the tiny 
grain case. When large grains and tiny grains coexist, the 
effect of large grains becomes negligible if xpah > 10~^. 

Although tiny grains strongly enhance PPD accre- 
tion in the outer disk, the situation at the inner disk 
is st ill similar to the case with 0.1 /^m grains (jBail 
l2011t ). with predicted accretion rate much less than 
the typical value of 10~^Mrr) yr~^ as inferred from ob- 
servations (jHartmann et al .[119981 ). The recently pro- 
posed far ultraviolet (FUV) ionization scenario does 
not provide large accretion rate in the inner disk ei- 
ther due to the small penetration de pth of FUV pho- 
tons (jPerez-Becker fc Chianej l2011bl ). Therefore, ei- 
ther additional ionization sources has to make the MRI 
more efficient than the X-rays and FUV photons (candi- 
date may include energ etic protons from the protostars 
(jTurner fc Drakell2009[ )). or additional mechanism such 



as magnetized wind (jSalmeron et al.l l2007| ) operates to 
provide the angular momentum transport in the inner 
disk. 

In the case of tra nsitional disks characterized by in- 
ner h oles or gaps () Gal vet et al.l 120021 : lEspaillat et al.l 
|2007() . we note that the observationally inferred outer 
boundary of the holes o r gaps is typically at a few 
tens of AU (|Hughes et al.l 120091 : iKim et al.l[2009( ). which 
is generally greater than rtrans in our models with 
PAHs. Therefore, the enhancement of accretion by tiny 
grains works effectively for transitional disks, and X- 
ray driven MRI with PAHs is able to feed the inner 
hole /gap at the (optimistic) rate of about 1O~^M0 
yr~^. This is sufficient to accoun t for the observed 
accretion rate in tra.nsitio nal disks (jNaiita et al.l 120071 : 
ISicilia-Aguilar etaLJ |2010[ ). with the accreting gas fed 
from the outer disk flowing through the h ole / gap possi- 
bly guided by multiple planets ([C hiang fc Murray- Clavl 
l2007tlPerez-Becker fc Chiandl2011aHZhu et al.ll20l'lir 
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APPENDIX 
MIDPLANE ACTIVATION AT TRANSITION RADIUS 

It has been shown in Figure |4] and discussed in Section 13.3! that near the transition radius rtrans where the disk 
midplane becomes active, the predicted accretion rate M increases sharply with radius, from well below the grain-free 
rate at r < rtrans to a factor of several higher at r > rtrans- Wc discuss the activation of the midplane by the MRI in 
this Appendix and explain the sharp dependence of M on disk radius near rtrans- 

Figure [5] shows the MRI permitted region in our fiducial model at 10 and 16 AU between which lies the transition 
radius. We see from this plot that much stronger magnetic field is permitted in the disk midplane once it is activated. 
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which is the cause of the big jump in M . 

Looking more closely into the sharp transition, we find it is related to the dependence of the AD Elsasser number 
Am on the magnetic field strength. We have shown in Figure [1] that rjA depends on the magnetic field quadratically 
in weak (/3e < 1) and strong (/3i > 1) field regimes as commonly considered, but at the intermediate regime with 
/3j < 1 < /3e, r]A behaves similarly as rjo and does not depend on the magnetic field strength. Consequently, the 
contours of Ohmic and AD Elsasser numbers A and Am arc nearly parallel to each other, and both shapes are close to 
(near the midplane) or slightly steeper (near the surface) than the contours of constant f3, as can be seen in Figure [HI 
In the mean time, an increase of ionization level (as one moves to outer radius) brings down the Am and A contours 
relative to the /? contour in a wide range of disk heights. Then, according to our criteria for the MRI to operate (A > 1 
and /3 > f3jnin{Am)), once the ionization level in the disk reaches some critical degree (i.e., left panel of Figure IH]), 
a further increase in ionization would suddenly make both criteria satisfied simultaneously across the midplane (i.e., 
right panel of Figure [6|), leading to the sharp transition. 
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